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We investigate dynamical symmetry enlargement in the half-filled SU(4) Ifubbard chain using 
non-perturbative renormalization group and Quantum Monte Carlo techniques. A spectral gap is 
shown to open for arbitrary Coulombic repulsion U . At weak coupling, f/ ^ 3t, a SO(8) symmetry 
between charge and spin-orbital excitations is found to be dynamically enlarged at low energy. At 
strong coupling, U > 6t, the charge degrees of freedom dynamically decouple and the resulting 
effective theory in the spin-orbital sector is that of the SO (6) antiferromagnetic Heisenberg model. 
Both regimes exhibit spin-Peierls order. However, although spin-orbital excitations are incoherent 
in the S0(6) regime they are coherent in the SO(8) one. The cross-over between these regimes is 
discussed. 

PACS No: 71.10.Fd, 71.10.Pm, 02.70.Ss, 71.30.-fh 



In strongly correlated electronic systems the presence 
of additional dynamical degrees of freedom out of the 
usual spin and charge ones is expected to play an im- 
portant role in a number of complex systems. This is 
the case, for example, of some d-electron systems [1], 
Ceo-based materials [2] and, also, various ladders-type 
compounds [3], for which low-energy excitations can- 
not be constructed from a single effective orbital per 
site (one band). An important question that arises 
is to know whether or not there exist generic fea- 
tures associated with multi-orbital effects. Such a ques- 
tion is non-trivial since it is known that the lack of 
symmetry in multi-orbital problems (beyond the usual 
SU(2) spin-invariance) is responsible for the presence of 
many independent couplings and, therefore, a wide range 
of problem-dependent physical behaviors could be ex- 
pected. However, at sufficiently low energy, it may hap- 
pen that the effective symmetry is increased, thus consid- 
erably simplifying the description of the problem. This 
is of course what happens in a critical (gapless) model. 
In the more general case where a spectral gap is present, 
the possibility of such a Dynamical Symmetry Enlarge- 
ment (DSE) at low energy is clearly non trivial. Recently, 
Lin, Balents, and Fisher [4] have emphasized that DSE is 
likely to be a generic tendency of the perturbative (one- 
loop) Renormalization-Group (RG) flow in their study 
of the half-filled two-leg Hubbard ladder. However, since 
in a gapped system DSE is a strong coupling effect one 
may thus question the reliability of perturbation theory 
[5]. Clearly, in view of the importance that such a DSE 
phenomenon might have in our understanding of com- 
plex systems, a non-perturbative investigation is called 
for and it is the purpose of this Letter to present such a 
study. 

In the following, we investigate the DSE phenomenon 



using non-perturbative RG and Quantum Monte Carlo 
(QMC) simulations for the simplest one-dimensional half- 
filled two-band Hubbard model where spin and orbital 
degrees of freedom play a symmetrical role. The corre- 
sponding SU(4) Hubbard model reads as follows: 



U 



(1) 



where c\ creates an electron with spin a = (t, i) and 
orbital index a = (1,2) at the ith site, and n^ aa ~ 
cl^^Ci^aa- The total symmetry group of (1) is U(4) — 
U(l)chargo<8) SU(4)spin-orbitai, it is the maximal symme- 
try allowed for a two-band Hubbard model. A simple 
one- loop perturbative analysis [4] would predict that, at 
half- filling, a SO (8) symmetry between charge and spin- 
orbital degrees of freedom is likely to be dynamically 
enlarged at low energy. Such a DSE pattern, U(4) 
SO (8), is highly non-trivial since one naturally expects 
the charge degrees of freedom to decouple at sufficiently 
large U. Indeed, in the limit U ^ t the Hamiltonian 
(1) reduces, at half-filling, to an antiferromagetic (AF) 
Heisenberg model (where the spin operators act on the 
six-dimensional antisymmetric representation of SU(4)). 
It is precisely the interplay between the small-?/ predicted 
SO (8) regime and the large-C/ charge-decoupled Heisen- 
berg limit which is considered here. 

The low-energy effective field theory associated with 
(1) is obtained, as usual, by performing the continuum 
limit and is most suitably expressed in terms of the exci- 
tations that are related to the symmetry of the problem. 
The U(l)chargc charge sector is described, in a standard 
way, by a single bosonic field $c and its dual field 9c. 
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There are many equivalent ways to describe the spin- 
orbital excitations in the SU(4)spin-orbitai sector and it 
is most convenient to represent them by six real (Majo- 
rana) fermions a = (1,...,6) [6]. We find that the 
low-energy effective Hamiltonian associated with (1) is 
given by: 



^= 2" 



-il%iV)+T^9s\Y.''"\ -2i5«cCos(v^$,)^K^ (2) 

\a=l / a=\ 

where = gs = -gsc = -U/2Tr, and Vg = vp+Qs, 

vp = 2t being the Fermi velocity. In Eq. (2), the Lut- 

tinger exponent Kc 



= 1/ y/l + 2gc/vF and the charge 
velocity Vc = vpy/^ + Ig^jvp depend on the charge cou- 
pling gc = "iUj^-K. The low-energy effective field the- 
ory (2) describes the interaction between a SO (6) Gross- 
Neveu (GN) model, associated with spin-orbital degrees 
of freedom, and a Luttinger liquid Hamiltonian in the 
charge sector. The interaction term, with coupling con- 
stant gsci is an umklapp contribution that comes from 
the 4ki^ part of the Hamiltonian density and is only 
present at half-filling fci? = 7r/2. In sharp contrast with 
the half-filled SU(2) Hubbard model and the SU(4) case 
at quarter- filling [7], there is no spin-charge separation 
at low-energy at half-filling. Spin-orbital and charge de- 
grees of freedom remain strongly coupled through the 
4ki? umklapp process. At this point it is worth stress- 
ing that there exists a higher-order umklapp term (8kf 
process) Vc = y cos(\/l67r$c) which depends only on 
the charge degrees of freedom. Although this operator, 
with scaling dimension A = ^K^., is strongly irrelevant 
at small [/, it may become relevant at sufficiently large 
U . As we shall see, this contribution is at the heart of 
the physics of the SU(4) Hubbard model in the large U 
limit. 

A simple one-loop RG calculation reveals that the cou- 
plings Qa = {gc,gs,gsc) flow at strong coupling. In par- 
ticular, as gc blows up at low energy Kc inevitably de- 
creases until Vc becomes relevant. Thus, the nature of 
the low-energy physics depends on the balance between 
the two umklapp operators with very different proper- 
tics. Clearly non-pcrturbative methods are called for. 
In this respect, Gerganov et al. [8] have provided an 
RG framework which allows to compute the RG /? func- 
tion to all order in perturbation theory for a large class 
of one-dimensional models with current-current interac- 
tions. We have applied their formalism to the Hamilto- 
nian (2) and obtained the resummed function. The 
detailed analysis of the non-pcrturbativc RG flow will 
be discussed elsewhere [9] and we shall here present our 
main result. Neglecting velocity anisotropy, we get: 
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where ija = dga/dt, t being the RG "time" and ga 
ga/vp- In absence of the umklapp contribution Vc, we 
find that the RG flow crucially depends on gc as follows. 
In the weak-coupling regime, at small enough U /t such 
that 5c < 2, all the couplings converge to the same value, 
9a{t*) = 2, at some finite RG "time" t*. On the other 
hand, when gc > 2, one enters a regime where perturba- 
tion theory is meaningless. 

Weak Coupling Regime. When gc < 2, much can be 
said on the low-energy physical properties of the model 
(1). Indeed, integrating the fiow up to t* , one finds that 
the Hamiltonian (2) at that scale reduces to the 80(8) 
GN model: 
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where we have refermionized the charge degrees of free- 
dom in terms of two real fermions ^'''^: {Sj + ii^)R(L) ~ 
exp (±z\/47r^ci?(L))- The equivalence at low energy 
between (1) and (4) is a manifestation of the DSE 
U(l)chargo® SU(4),pin_orbitai ^ S0(8). This S0(8) en- 
larged synunetry which has been first predicted using a 
1-loop RG calculation in [4] , is shown here to hold beyond 
perturbation theory provided gc <2. For higher values of 
gc, the higher-umklapp term Vc plays a prominent role at 
low-energy and, as we shall see. is responsible of the dy- 
namical decoupling of the charge degree of freedom. One 
of the main interest of the emergence of this SO (8) sym- 
metry stems from the fact that the model (4) is integrable 
and a large amount of information can be extracted from 
the exact solution [4,10]. The low-lying spectrum of the 
SO (8) GN model (4) is fully gapped and consists of three 
distinct octets with the same mass m ~ The 
fundamental fermion octet, associated with the Majo- 
rana fermions of Eq. (4), is made of two charged 
±2e spin-orbital singlets, called cooperons, and six spin- 
orbital excitations which transform according to the self- 
conjugate representation of SU(4) with dimension 6. The 
remaining two octets arc of kinks type. In particular, the 
excitations of the SU(4) Hubbard model (1), carrying the 
quantum numbers of the lattice fermions Ci^aa, are rep- 
resented by eight of these kinks. In addition, there are 
28 bosonic states organized as a rank-2 S0(8) antisym- 
metric tensor and a singlet, all of mass -\/3m which can 
be viewed as bound states of the fundamental fermions 
or of the kinks states. The massive phase corresponding 
to the S0(8) GN model (4) is a spin-Peierls (SP) phase 
as it can be readily shown by considering the order pa- 
rameter Osp = J2i,aa i~^y4,aa^t+i,aa which has a non- 
zero expectation value (Osp) ^ 0. The ground state of 
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(4) is thus doubly degenerate and spontaneously breaks 
the lattice translation symmetry. The striking feature 
is that, despite of this dimerization, both electronic and 
spin-orbital excitations are coherent, i.e. they contribute 
to sharp peaks in various spectral functions. This result 
stems from the existence in the exact spectrum of (4) of 
states that have the same quantum numbers as the elec- 
tron and spin-orbital operators cl^^Ci^bfS- In particular, 
the dynamical structure spin factor of the system displays 
a sharp peak at energy lo = \f2>m corresponding to an ex- 
citation of one of the bosonic states of the SO (8) theory. 
In this respect, the half- filled SU(4) Hubbard model is 
predicted to be a fully coherent gapped dimerized liquid 
at weak-coupling. 
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FIG. 1. Onc-particlo gap Ai, spin-orbital gap As, and 
cooperon gap Ac as a function of U/t; inset: A^ gap as func- 
tion of U/t. 

At this point it is worth discussing the stability of this 
SO (8) phase. Although the RG equations (3) are non- 
perturbative in nature it remains to investigate the ef- 
fect of neglected symmetry-breaking operators such as 
the higher-umklapp term Vc and chiral interactions that 
account for velocities anisotropy. For small symmetry- 
breaking terms, the SO (8) multiplets will be adiabatically 
deformed and split into U(l)chargc x SU(4)spin-orbitai 
multiplets: the SO (8) symmetry is only realized approxi- 
mately at weak enough coupling. At small U/t the split- 
tings are exponentially small but we expect perturbation 
theory to break down as U increases even when Qc < 2. 
The reason stems from the neglected umklapp operator 
Vc which becomes relevant before one reaches the SO (8) 
symmetry restoration point as A < 2 when gc > 3/2. We 
thus expect the SO (8) regime to hold approximately up 
to some critical value Uc of which a very naive estimate 
can be obtained using the bare value of gc- Uc'^ 2TTt. 

In order to check our theoretical predictions we have 
performed extensive T — QMC simulations of the 
SU(4) Hubbard model (1) at half filling for a wide range 
of U /t. Following the work done in Ref. [7] in the quarter- 
filled case, we have computed all gaps associated with 
the S0(8) tower of states. We discuss here our results for 



three of them: Ai which is the gap to the one-particle ex- 
citation c]^„^, Ag which is the spin-orbital gap associated 

with the excitations c-^^Cj^b^ and finally the cooperon 
gap Ac which is the gap to a spin-orbital singlet state of 
charge 2e. The latter excitation is a striking feature of the 
SO (8) spectrum and is not simply related to electronic 
excitations on the lattice. For example, the cooperon 
comes into pairs from the charge 4e excitation ^aaclaa- 
We have computed the cooperon gap Ac as half the gap of 
this state. The exact spectrum of (4) imposes the highly 
non-trivial predictions for the ratios: (Ai/Ac)so(8) = 1 
and (A,5/Ac)so(8) = \/3- Strong deviations from these 
theoretical predictions will be a signature of the failure of 
the increased S0(8) symmetry. We show in Fig. (1) our 
results for Ai(C/), As{U) and Ac{U) for values of U/t 
ranging from 0.5 to 20. The extrapolation to the ther- 
modynamical limit has been performed using lattice sizes 
L = 8, 16, 32, 48, 64 and the errors on the gaps range from 
10^^ at small U/t to 10~^ at large U/t. Two asymptotic 
regimes are identified. A small U/t regime and a large 
U/t regime where spin-orbital and charge degrees of free- 
dom clearly separate. Both regimes are most easily seen 
on the spin gap As(f7) behavior (inset of Fig. 1) which 
increases until it reaches a maximum around U/t ^ 6 and 
then decreases smoothly to zero asU/t oo. Clearly the 
SO (8) regime is expected to show off at small U/t. 
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FIG. 2. Gap ratio Ai/Ac as a function of U/t. 

In Fig. 2 we plot the ratio {Ai/ Ac){U/t). Despite our 
high quality QMC datas, we have not been able to re- 
solve the ratios at very small U's where the gaps are 
exponentially small. However, one observes a clear sat- 
uration of the ratio at the SO (8) value as U decreases 
below U ~ 3.5t. Other gap ratios, not presented here 
[9], show also a S0(8) saturation in the regime U < 3.5*. 
These results strongly support the existence of a SO (8) 
DSE in a regime where the gaps are not infinitesimally 
small (As ~ O.li - 0.2t). Above U ~ 3.5t, the ratio 
shows a departure from its SO (8) value. Though such a 
behavior may be attributed to level splitting due to sym- 
metry breaking operators at small U/t, this is certainly 
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not the case above U ^ 6t where Ai/Ac saturates at the 
value 1/2. It is difficuh from our results to give a pre- 
cise value to Uc above which the S0(8) regime is lost but 
we can give an estimate 3t < Uc < 6t. Notice that the 
upper value is in agreement with our rough estimate of 
Uc = 2nt based on a scaling argument. We shall see now 
that the physics at large U/t is of a very different nature. 

Strong Coupling Regime. When U/t ^ 1 there is 
a clear separation between spin-orbital and charge de- 
grees since Ac Ag (see Fig.l). The umklapp term Vc, 
which depends only on the charge degrees of freedom, 
becomes now much more relevant than the Akp coupling 
and charge fluctuations are strongly suppressed by this 
process. Integrating out the charge degrees of freedom, 
the low-energy effective Hamiltonian in the spin-orbital 
sector reduces to a massive S0(6) GN model: 

6 6 
a— 1 a=l 

+ ^«(^)(E^") , (5) 

where M > and Gs{U) is a negative effective coupling 
at large U/t. The Hamiltonian (5) describes six massive 
Majorana fermions with a weak repulsion. One can show, 
using Eq. (5), that (Osp) 7^ so that the ground state 
is still in a SP phase. Neglecting charge fluctuations, 
this dimerized phase, with broken translational symme- 
try, can be simply understood as a set of nearest-neighbor 
SU(4) ^ SO (6) spin-orbital singlet bonds. There is thus 
a continuity between weak and strong coupling with re- 
spect to the nature of the ground state. However, there 
is a striking difference between the SO (8) regime and this 
strong coupling phase, called SO (6) regime, at the level of 
the coherence of excitations. The excitation spectrum of 
the model (5) for < consists of massive fermions 
which are the SU(4) dimerization kinks, and their mul- 
tiparticle excitations. In particular, there are no bound 
states so that the spin-orbital dynamical structure fac- 
tor exhibits a two-particle continuum: the spin-orbital 
excitations are now incoherent. Apart from these neu- 
tral excitations, there are massive modes corresponding 
to solitons in $c with charge q = ±e coupled with zero 
modes of the Majorana fermions of Eq. (5). These 
excitations have a larger gap and carry the same quan- 
tum numbers as the kinks of the S0(8) spectrum [9]. 
The cooperon is no longer a stable excitation in the large 
U/t limit but becomes instead a diffusive state made of 
these two kinks. One thus expects that the gap ratio 
Ai/Ac saturates at 1/2 in the SO (6) regime in full agree- 
ment with the numerical results of Fig. (2). The physics 
of the strong-coupling regime can also be investigated 
by a complementary approach which consists to map di- 
rectly, in the large U limit, the SU(4) Hubbard model (1) 
onto a SU(4) AF Heisenberg chain by a standard pertur- 
bation theory in t/U [11]: HeS = J^i^i ■ ^i+i, with 



J = At^ /U and Sf' are SU(4) spin which belongs to the 
six-dimensional representation of SU(4). This S0(6) AF 
Heisenberg chain is not integrable and has been studied 
by means of the density matrix RG approach [11]. In 
full agreement with our results, this model belongs to a 
SU(4) dimerized phase. Using the numerical results of 
Ref. [11], we find that our QMC results for the spin gap 
As{U/t) follow a S0(6) Heisenberg regime for U > 8t. 
In this respect, we deduce that the low-energy physics of 
the S0(6) AF Heisenberg chain is described by the six 
almost free massive Majorana fermions (5). 

Cross-Over Regime. Both S0(8) and S0(6) regimes 
differ by the coherent nature of spin-orbital excitations 
and the existence of an elementary charge 2e cooperon 
excitation. In the simplest hypothesis, the cross-over be- 
tween these two regimes can be understood as a change 
of sign of the coupling Gg as a function of U. A mean- 
field analysis of the low-energy effective theory together 
with our numerical result predicts that such a cross-over 
occurs at [/ ~ 4.5t [9]. When Gs{U) > 0, the Majorana 
fermions of Eq. (5) experience an attractive interaction 
and neutral bound-state in the adjoint representation of 
SU(4) are formed. The latter excitation is adiabatically 
connected to one of the bosonic states of the SO (8) spec- 
trum which is responsible of the sharp peak in the dynam- 
ical structure factor in the SO (8) regime. It is thus very 
tempting to conclude, within this simple scenario, that 
the SO (8) regime approximately extends up to Uq ^ 4. St 
above which one enters the Heisenberg S0(6) regime. 
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